Quantum Nyquist Temperature Fluctuations 
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We consider the temperature fluctuations of a small object. Classical fluctuations of the tem- 
perature have been considered for a long time. Using the Nyquist approach, we show that the 
temperature of an object fluctuates when in a thermal contact with a reservoir. For large tempera- 
tures or large specific heat of the object C v , we recover standard results of classical thermodynamic 



fluctuations (AT ) = 



Upon decreasing the size of the object, we argue, one necessarily 




reaches the quantum regime that we call quantum temperature fluctuations. At temperatures be- 
low T* ~ fi/fesT, where r is the thermal relaxation time of the system, the fluctuations change the 
character and become quantum. For a nano-scale metallic particle in a good thermal contact with 
a reservoir, T* can be on a scale of a few Kelvin. 

PACS numbers: 05.40.-a, 07.20.Dt 



Classical thermodynamic fluctuations have been stud- 
ied for more than two centuries. The Gibbs canonical dis- 
tribution function, P(a)da oc exp[—E(a)/kBT]da, is one 
of the most fundamental concept in statistical physics. 
All thermodynamic variables can be obtained from this 
distribution function. According to the initial assertion 
of Gibbs, the temperature of a canonical ensemble is con- 
stant and thus does not fluctuate. Therefore, the tem- 
perature fluctuation cannot be generically represented by 
the above distribution function. Instead it is derived 
from energy fluctuations (SE). Alternatively, the von 
Laue approach [EJ to fluctuating system thermodynam- 
ics via the minimal work can also lead to the temperature 
fluctuation (ST) . In recent years, there has been increas- 
ing interest in the nano-scale problems such as the glass 
transition 0, nucleation ||, and protein folding Of 
device importance, the mechanical resonators are being 
pushed to the nanometer scale J|, |(| . For these nano-scale 
systems, the temperature fluctuation can be large. Re- 
cently, it has been shown J?J that the von Laue approach 
gives much more reasonable results for the temperature 
fluctuation in a confined geometry than the treatment 
with the Gibbs distribution function. To the best of our 
knowledge, the existing study has been limited to the 
classical regime. We know that any classical variable, 
say coordinate, force, etc., has its corresponding stan- 
dard quantum limit where quantum fluctuations domi- 
nate. Similarly, we expect temperature T will have its 
quantum limit. Here we argue that when the tempera- 
ture is below T* ~ H/t, where r is the thermal relaxation 
time of the nano-scale particle, a quantum temperature 
fluctuation regime emerges. 

Consider an experiment in which we are going to mea- 
sure the temperature fluctuation of an ensemble of in- 
creasingly small objects. Without loss of generality con- 
sider a set of quantum dots, as shown in Fig. [I]. Assume 
these dots to be similar in the number of contained par- 
ticles, size, etc. In addition, each dot has discrete levels, 
which arc filled by a sufficient number (N 3> 1) fermions 
(e.g., electrons) or bosons (e.g., 4 He atoms). All these 
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FIG. 1: The ensemble of the quantum dots. The thermal 
relaxation time within each dot is tj.. Each dot is thermally 
coupled to the reservoir with temperature T. Note that each 
dot contains large enough number of particles that the ther- 
modynamics consideration is still applicable. In particular 
the notion of a quasi-equilibrium distribution is applicable to 
each dot. Consider a thermal fluctuation that produces the 
temperature fluctuation in the ith dot, T + STi. After the 
equilibration time Ti , the dot relaxes to the quasi-equilibrium 
state with temperature T. Repeating this measurement on 
a set of dots, we find the distribution of temperature in an 
ensemble of dots. The width of this distribution of T will 
determine the fluctuation of temperature ST. 



dots are in contact with a substrate (large plate) which 
plays the role of a thermal reservoir. The reservoir is 
kept at a certain temperature T. The thermal contact 
between the dots and the reservoir will cause the ther- 
mal fluctuations in the dots. As a result, the heat flows 
to/from the reservoir. The relaxation time for the ther- 
mal process between the dots and the reservoir is r. 

Let us first briefly recall the derivation of the tem- 
perature fluctuation from the von Laue approach [jO. 
The propability w for a fluctuation is proportional to 
c?cp(—R m i n /kBT), where R m i n is the minimum work 
needed to fulfill reversibly the given change in the ther- 
modynamic quantities in the quantum dot and ks is 
Boltzmann's constant. For simplicity, we assume the vol- 



ume of each dot be fixed so that R n 



AE - TAS, 



where AE and AS are respectively the changes in the 
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energy and entropy. Therefore we have: 

. AE - TAS\ 
w = AIq cxp[- 



(1) 



where A/"o is the normalization constant. For small fluc- 
tuations, by expanding AE to the second order in AS, 
and noticing that AS = (C V /T)AT with C v the specific 
heat, it is found: 



w = Afi exp[- 



C V AT 2 



(2) 



2k B T 1 1 

where TVi is a new renormalization constant. It follows 
immediately that the average square fluctuation of tem- 
perature at a constant volume is: 



(AT 2 ) = T 2 /C v 



(3) 



This is the result from the standard classical theory of 
thermodynamic fluctuations. However, it has to be kept 
in mind that for the starting equation (Q) to be valid, 
the temperature has to be much bigger than the thermal 
relaxation rate, i.e.: 



k B T > H/t 



(4) 



which also means when the temperature is too low or r 
is too small, the fluctuations can no longer be treated 
classically. 

Here we propose the use of the Nyquist approach M to 
treat the temperature fluctuation. First consider a gen- 
eralized coordinate x and its relaxation % = — X(x — xq), 
where Xq is its equilibrium value. The external force F, 
conjugated to the coordinate x, determines the equilib- 
rium xq. Now, in case there are fluctuations in the exter- 
nal force F, equilibrium value xq = xq + Axq fluctuates 
around its steady state position Xq by Axq = ^\ So AF. 
For the equation of motion we get 



dx 

— = —X[x 
dt K 



x 



(5) 



Assume now temperature T = x to play the role of a 
generalized coordinate and the entropy S = F the role 
of the generalized fluctuating force. The relaxation pro- 
cess of the temperature can be described by a linearized 
macroscopic "equation of motion" : 



r = -A(Ar-g b A S ), 



dt 



(6) 



where A = 1/r, and AT = T — To, and is the deviation of 
the equilibrium temperature T as a result of the fluctuat- 
ing force AS. Equation (Q) is valid for the positive time 
and can be extended to the negative time by changing 
sign of the derivative. Performing the Fourier transform 
for AT 



AT(t) = 



1 

2^ 



dioAT^e 



(7) 



and similarly for AS, we can arrive at 
ATL - a(w)A5 aJ , 



(8) 



where the response function or generalized susceptibility 
reads: 



a(uj) 



XT 



C v (—iuj + A) 



(9) 



where we have used j^\t = Using the fluctuation- 
dissipation theorem as developed by Callen and Wel- 
ton, which relates the fluctuation of a thermodynamic 
quantity to the imaginary part of the susceptibility 
a(ui) jl], [|, H^|, it immediately follows: 



(AT 2 ) U = hcoth(haj/2k B T)a"(uj) 
where the imaginary part of a(to): 

XT u 

a (lu) 



C v u 2 + X 2 



(10) 



(11) 



For the average quadratic fluctuation of T, it can be 
found: 



(AT 2 ) = 



HXT 
2ttC v J_ 



du>- 



uj 2 + A ; 



■coth(/kj/2fc B T) . (12) 



The integral on the right side of Eq. (|12| ) depends on 
the ratio of ksT/hX. When ksT/HX ^> 1, by expanding 
coth(hw/2k B T) = 2k B T jhio + hio /(Sk B T for hw/2k B T <C 
1, we have 



(AT 2 ) 



k B T 2 
C (. 



1 



HX hljj r 

■ In ■ 



irk B T k B T 



(13) 



where we have introduced an upper band cutoff oj c ~ 1 /t± 
on the order of the relevant bandwidth since at the high 
frequency the integral J Q c dw ji^^i = hi(Lu c /rj) is loga- 
rithmically divergent as w c — > oo. One can recognize 
immediately that Eq. ( |l3| ) is the classical limit of the 
temperature fluctuations as derived above from the von 
Laue approach. In the opposite limit of low tempera- 
tures, hX ^> k B T, one finds: 



(AT ) = 



? . hXT oj c 
In — 



ttC v X 



(14) 



Therefore, we find that at low temperatures the tem- 
perature fluctuations would acquire a distinctly quan- 
tum character with H/t entering into the magnitude of 
(AT 2 ). From ergodicity assumption, it follows that the 
time averaged temperature of each particular dot is equal 
to the average value of T. Any fluctuation, described by 
Eq. (pT|) , happen on a characteristic time scale r. 

The high temperature expansion in Eq. ([TJ]) has al- 
ready indicated the crossover temperature 



HX LU r 

T ln-T 

K B TT X 



(15) 
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FIG. 2: Temperature dependence of the fluctuations ST. As 
r>T*, the classical regime is recovered. As the temperature 
gets smaller or the subsystem is getting smaller, a quantum 
fluctuation regime comes into play. Note that the specific 
heat C v (~ T for fermions and ~ T 3 for phonons at low 
temperature) is temperature dependent. 



at which there is a change of the regime from the clas- 
sical to quantum fluctuations. Physically, T* sa h/ksT 
corresponds to the uncertainty in energy associated with 
the relaxation process in the subsystem. The reservoir 
is attached to a subsystem via a thermal contact that 
has its own bandwidth H/t and any temperature fluctu- 
ation will relax on the scale of r. Once T -c T*, the 
intrinsic bandwidth of the contact rather than the tem- 
perature will dominate the Gaussian fluctuations. As 
stressed in Ref. jy , fluctuations cannot be treated classi- 
cally if a fluctuating quantity is changing too rapidly or 
if the subsystems are too small. 

Our derivation, which is essentially identical to the 
classical to quantum Nyquist crossover in the standard 
Nyquist theory, provides the description of a quantum 
regime of the temperature fluctuations. There are few 
possible limitations of this description. One is that our 
approach in not applicable to the systems outside of ther- 
mal equilibrium, say hot electrons and cold phonon bath, 
as sometimes is the case. Another restriction is that the 
typical energy level spacing 5 ~ ^/L d , d is the dimen- 
sionality of the dots, should be small compared to T. If 
the temperature is much smaller than the level spacing, 
there will be no thermally excited states and the notion of 
the thermodynamic equilibrium is not applicable to this 
system. We should also point out that for small particles 
there are mesoscopic corrections to the total energy of 
particle, that can be related to the temperature fluctua- 
tions. These mesoscopic fluctuations have a typical tem- 
perature scale given by Thouless energy ksTrh = D/L 2 , 
L - is the typical size of the particle, D- is diffusion co- 
efficient, and occur in addition to the fluctuations we 
consider here. We assume here that T* > Txh- 



Experiments on the temperature-dependent fluctua- 
tions of magnetization of small paramagnct were per- 
formed by Chui et al. [jll| . Experimentally observed spec- 
tral density (AT 2 ) W was shown to be of the form given by 
the high temperature expansion of Eq. (|lOj). The ther- 
mal relaxation time was r ~ 1 second for the considered 
size of the paramagnet (about I cm 3 ). The total tem- 
perature fluctuation (AT 2 ) = T 2 /C v was also claimed 
as a result of the integration over Eq. @. Therefore 
we can regard this result as an experimental evidence for 
the classical temperature fluctuation in the canonical en- 
semble according to Eq. (|l0|). The obvious next step is 
to extend these measurements to the samples of much 
smaller sizes, down to 1 pm in its linear size and study 
the temperature dependence of the fluctuation spectrum 
at low temperatures. 

Now let us estimate the crossover temperature T* for a 
metallic dot and for a 4 He droplet. The low temperature 
limit T <C hX/ks implies that the relaxation time of the 
thermal object has to be short enough. Since the thermal 
relaxation time r = C v Rt, where Rt is the thermal resis- 
tance of the contact between the object and the thermal 
reservoir. The easiest way to achieve a short thermal re- 
laxation time is to work with the smaller objects (smaller 
C v ) with good enough thermal contact with the reservoir. 

(i) For the metallic system, we consider a nanometer 
mechanical resonator, which is a cylindrical gold (Au) 
rod of 1 yitm in length L and 15 nm in radius r. The mass 
of this Au rod is tuau = PAuLA — 1.36 x 10 -14 g, where 
the mass density of Au pAu — 19.3g/cm 3 , and A = nr 2 
is the cross-section area. The mole mass of Au is 196.97 
g/mole, we thus find the specific heat constant for Au 
per gram is -f A u = 0.73mJ/mol • K 2 = 3.7 x I(T 6 J/gK 2 . 
At T = 0.1 K, C^ u = rriAujAu = 4.76 x Ifr 21 J/K. The 
thermal conductivity of Au at T = 0.1 K is about Kt = 1 
Watt/K m. The thermal resistance is Rt = -^p^ — 
1.4 x 10 9 K • Second/J. One then finds r ~ 6.7psec, and 
T* ~ 1 K, respectively, which is now experimentally ac- 
cessible. In practice, the thermal impedance mismatch at 
the interface between the nano-scale subsystem and the 
reservoir would lead to a much lower conductance and a 
longer r. To estimate the role of the "bottleneck", one 
would need a specific model. However, from the above 
estimate, one gets an impression that, regardless of the 
experimental constraints, the quantum temperature fluc- 
tuation below T* is observable. 

(ii) For the case of a small bosonic system, we consider 
a droplet of 4 Hc, which is enclosed in a metallic con- 
tainer such as lead. The size of the droplet is taken to 
be 0.1 pm. Below the super/fluid transition temperature, 
the specific heat of 4 He is dominated by phonons, which 
follows the power law as 0.02 x T 3 J/gK 4 Q. For the 
above given size of the droplet and at T ~ 0.3 K, the spe- 
cific heat C v = 2.7 x 10~ 18 J/K. At these temperatures, 
the thermal resistance is dominated by a surface resis- 
tance due to the contact of the droplet with the metal. 
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From Fig. 8.6 in Rcf. (12|, we estimate the thermal resis- 
tivity Rt ~ 2Kcm 2 /Watt between the 4 Hc droplet and 
the metal surface. One then obtains the relaxation times 
t ~ 6.9 x 10 -8 second, which corresponds to T* ~ 10 _3 K, 
which is small for the given size of the droplet. 

Experimentally the proposed crossover to quantum 
regime can be seen as a change in temperature depen- 
dence of noise of some observable. The choice depends 
on a specifics of the experiment obviously, e.g for an os- 
cillating clamped beam it can be a noise of the me- 
chanical oscillatior. In the case of magnetization noise 
PI , one would desire to measure noise in the SQUID at 
relevant frequencies 1/r. 

In summary, we have used the Nyquist approach to 
study the temperature fluctuation of an object in ther- 
mal contact with a reservoir. It is shown for the first 
time that when at temperatures below a characteristic 
value T* ~ h/ksT, the temperature fluctuation would 
acquire a distinctly quantum character. For a nano-scale 
particle, T* is on the order of a few Kelvin. In light of 
recent advances in nano-technology, the quantum fluc- 
tuation regime should be experimentally accessible and 
might be relevant for the experiments on nanoscale sys- 
tems. 
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